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Is the new state of matter formed in relativistic heavy ion collisions the 
deconfined quark-gluon plasma? We survey the status of several strange 
hadron observables and discuss how these measurement help understand 
the dense hadronic matter. 
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1. Introduction 

Matter made of free quarks has been discussed in literature for more than 
30 years. For further theoretical details and historical developments we refer 
to the 'Theoretical Foundations' collection of QGP articles f]]. Today, many 
believe that a new state of matter made of free quarks has been formed in 
relativistic nuclear collisions. Is this 'really' the interacting plasma of QCD 
quanta, quarks and gluons (QGP), formed and present in a limited domain 
of space and time, where quarks and gluons are propagating constrained by 
external 'frozen vacuum', which abhors color? Before proceeding to more 
specific matters, we will make a few general remarks. 

We shall discuss in these notes the yields of a diversity of strange hadronic 
particles produced in relativistic heavy ion collisions, and consider also theo- 
retical models of strangeness production in QGP. In the analysis of (strange) 
hadron yields we use the statistical hadronization model (SHM) to connect 
properties of matter with yields of particles (or vice- versa). We note that 
experiments detect reliably just a small fraction of all hadrons of interest. 

* Lecture Notes for the Cracow School of Theoretical Physics, XLVIe Course, 2006. 
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In order to have global and detailed information about the properties of the 
dense matter source of these hadrons, an extrapolation must be made to 
cover production yields of all particles. 

The yield of measured particles are obtained from their spectra, which 
are extrapolated to allow the integration over the full phase space. This 
introduces considerable systematic uncertainty, but reduces greatly the de- 
pendence on the dynamics of fireball evolution at time of hadronization. In 
this work, we address quantitatively only the mj_-integrated total, or rapid- 
ity density dN/dy, particle yields, noting sometimes the spectral m± shape. 
Even in this reduced case, considering the complexity of hadronization, and 
the very rich experimental data set, considerable effort is required to com- 
plete the data analysis. For this reason, we helped develop a hadronization 
tool SHARE (Statistical HAdronization with REsonances), a suite of nu- 
meric analysis programs which is available in the public domain 2 . SHARE 
is used in all analysis presented here. SHARE encompasses all other statis- 
tical hadronization models (SHM), which arise as special cases in SHARE 
set of parameters. 

As the above implies, there are several different approaches one can 
take in the study of particle yields using SHM. In our opinion, SHM is a 
physics-motivated model extrapolating from a subset to all particle yields. 
Therefore, we take the most elaborate version of the SHM model to compare 
with experiment: the results we present are obtained allowing for the full 
chemical non-equilibrium in the analysis of hadron yield data. 

We present here the key strangeness QGP observables and look at the 
results with two questions in mind: has the experimental result been pre- 
dicted to be a QGP consequence? Is an alternative explanation of the data 
available today, which does not contradict the behavior of the data being 
explained? We will also discuss how understanding of strangeness produc- 
tion at SPS and RHIC helps to prepare for the LHC energy range, and the 
low energy RHIC run. 

We address in turn the A/p (next section), Strange (anti)baryon en- 
hancement at SPS and RHIC (section I3J), The K + /-7r + -horn (section |IJ). 
In this context, we expand on the observation of a change in the reaction 
mechanism, favoring baryons at energies at or above the K + /-7r + -peak. We 
then introduce the bulk observable, strangeness per entropy s/S in section 
13 and show how it helps evaluate the ratio of degrees of freedom present. 
We then show how s/S can be computed at RHIC and LHC and present 
results showing independence of the result from initial conditions other than 
entropy content dS/dy. In section |7| we discuss strange hadron resonances 
and their importance to the diagnosis of QGP. We close the discussion, in 
section IS] with a short outlook at the developments in the future. 
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2. Ratio of A with p at SPS 

The first proposed strangeness and hadron yield ratio signature of QGP 
has been the relative yield of A and p [3J|1]. After cancellation of combinato- 
rial and phase space factors, this ratio is determined by relative quark yields 
available at hadronization. If no QGP were formed, one could at best hope 
for chemical equilibrium yields in the hadron gas matter, but especially at 
relatively low reaction energies this requires a magic touch, the chemical 
equilibration of strangeness in hadron born reactions. 

For both particles considered we include, aside of directly produced A, 
p, the yields of more massive resonances which decay into these particles. 
SHARE accounts for this important effect. However we present first the 
more intuitive historical line of thought: we assume that the resonance 
contributions multiply both A and p by similar enhancement factors. This 
being the case, resonance effect can be ignored at first when considering the 
A/p ratio. 

In a baryon-rich QGP environment, the light antiquark u, d abundances 
are suppressed. This is easily understood as result of u, d annihilation on 
u, d excess present where baryons are present. In the statistical hadroniza- 
tion model, this effect is described by the baryochemical potential ^ = 
3(/i u + ^d)/2 = 3fj, q . However, the strange antiquark yield s is also sup- 
pressed by the strange quark mass. Integrating the particle phase space 
characterized by a production temperature Th (hadronization temperature) 
we find: 



qgp ****** 7 9 QGP2T /? 

OGP 

7s i-„Q GP „QGPw T , , 

= ^GP°- 9e ' W 

7g 

The last equality follows for the currently accepted value m s /T/, ~ 0.7. 

The relative yield originating in the hadron phase comprises, in place of 
strange quark mass suppression, the hadron phase space suppression factor: 



HG / \ 3/2 
ll_ ( ) e -(^-m p )/T fe (u^-^)/T f 



HG fq 

HG 



7^G \m p J 



= 2|_l. 3e -180MeV/T /e (^G_ M HG )/T/ (2) 
7<J 

In the range Tf ~ 160 ± 20 MeV, the HG yield is significantly smaller com- 
pared to the relative yield from QGP. Moreover, in HG the multiplicative 
factor 7s /7g, will be well below unity at low collision energy. As a side 
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remark, note that the analysis of the experimental data suggests that this 
chemical equilibration factor reaches unity for most central and most ener- 
getic SPS reactions which is evidence for a new production mechanism of 
all particles. 

An important and self-evident result is that, using hadron phase space, 
it is indeed impossible to ever obtain for A/p a value that exceeds unity. 
However, the situation is different when QGP hadronizes. Let us recall the 
equilibrium hadronization conditions: 

T h ^T f , (s + s)^^(s + s) R , S^^S*, M?-M? = ^-A*?, (3) 

here S is the entropy and both s and 5 can be understood to be rapidity den- 
sity. We abbreviate above Q for QGP and H for HG. Moreover, irrespective 
of question of phase transition, for the purpose of comparison of production 
yields one normally assumes that chemical potentials are measured and thus 
yields are compared at a given value obtained from experimental data: 

fJ-u- fJ-s = Const. . (4) 

In either case much of enhancement of A/p (and other strange antibrayon 
relative yield in QGP compared to HG) is due to the reduced threshold of 
strangeness in QGP compared to HG. 

This argument correctly describes the yields of dominant fractions of 
(strange) particles which carry most of strangeness. For the rarely produced 
particles which antibaryons are at low energy the situation is more involved. 
At the onset of global hadronization, /i? GP — 0. This is so since in QGP at 
first s — s = 0, and strangeness is not clustered in baryon number objects 
in plasma. Considering that the HG phase is asymmetric at a given baryon 

OOP 

number, there is a buildup of /v with time, reversing with time the 
asymmetric emission of hadrons. 

For conditions we consider, the preferentially emitted particles are car- 
riers of s-quarks. A (and (qs) kaons ) is emitted at first with preference. As 
time progresses this leads to an increase in fif G in plasma phase such that 
symmetric emission and later s-quark emission dominance ensues, assuring 
that strangeness conservation s — s = is achieved among all produced 
particles by the end of hadronization. However, it has also been speculated 
that the early asymmetric emission of particles with s > s could lead to the 
formation of a strangeness enriched residual matter jS] . Extensive searches 
for this effect failed to confirm this reaction model. What is perhaps instead 

OCP 

happening is that at some point fig is large enough such that emission 
of particles with s > s dominates. 

At the beginning of the process of hadronization, the production of A 
(and H, Q) is enhanced by /i? GP = 0. Once produced, particles do not disap- 

OOP 

pear, and thus QGP property /v = is imprinted on the enhanced yield 
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of rarely produced strange antibaryons, which are predominantly produced 
in the early time stage of QGP hadronization. Thus, 



^0.9A? GP , \$ GP = e^ ' T »\ (5) 

QGP 



OCP 

It is well understood that the value of quark fugacity Au at hadronization 
increases with decreasing reaction energy. This implies that with decreasing 
reaction energy the relative yield A/p increases, which at first sight is an 
unexpected result. 

Beyond this qualitative discussion, originally presented more than 25 
years ago, a quantitative prediction which includes all resonance decays is 
shown in Fig. 47 of our Acta Physica Polonica review article from 1996 0, 
which is reproduced at hadronization with figure caption, see figure ^ The 
energy scale E/B (fireball energy content per baryon) is somewhat unusual, 
yet it corresponds almost exactly to y / SNN, the center of momentum frame 
nucleon pair reaction energy. Namely, for stopping of energy being 50% of 
participant stopping, the fireball energy content per baryon E/B, seen in 
figure E would be just the CM energy per colliding nucleon pair. Study of 
energy and participant number content in the fireball matter formed at low 
energy SPS reactions suggest that 50% relative stopping is the right order 
of magnitude. 

The result seen in figure ^ is not surprising, following on our qualitative 
discussion, we see how the ratio A/p is increasing with decreasing reaction 
energy. It is of some interest to compare this result with the AGS and the 
SPS experimental results. The data analysis has been evolving. The results 
are shown in figure [2 based on compilation of data and theoretical results 
by the NA49 collaboration 7 . We see that the central rapidity ratio A/p is 
well above unity at all available reaction energies. With decreasing reaction 
energy, this ratio increases, just as it is expected in our above discussion. 
The fact that to lowest considered collision energies and smallest reaction 
systems at AGS this trend persists (within the quite large error bars) could 
suggest that QGP formation followed by slow and continuous hadronization 
is the reaction mechanism governing the production of strange antibaryons 
at relatively low AGS energies. 

In the bottom of figure [21 the bulk matter predictions are shown, based 
on global SHM fits to the experimental data carried out by different groups 
(including our own results), under differing scenarios, all theoretical results 
have been assembled by the NA49 group and we refer to this work for 
further discussion and references. Comparing data with bulk yields we 
see that the additional enhancement of A/p we attribute to dynamically 
evolving fi s is clearly most pronounced at lowest reaction energies. In fact, 
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Fig. 1. Strange antibaryon ratio A/p, as function of E/B in a QGP-fireball for 
7 S = 1; solid lines are for full phase space coverage, short dashed line for particles 
with p± > 1 GeV and long dashed line for particles with m± > 1.7 GeV (reproduced 
from Acta Physic Polonica 



for the highest reaction energies, where sudden breakup of the QGP fireball 
is assured, the discrepancy between bulk hadronization and experimental 
data could perhaps be accounted for by conventional means. Specifically, 
the experimental A/p ratio needs to be reduced to correct for the included 
weak decays H — > A which are not included in the SHM models. Moreover, 
the thermal rapidity distribution of A is narrower than that of p, which 
enhances the central rapidity ratio compared to the 47r-yields shown in the 
theoretical part of the figure. On the other hand, the trend of the data at 
low energies is quite clear, despite the large error bars. 

We conclude: A/p > 1 ratio, rising with decreasing reaction energy, 
has been predicted to arise in hadronization of QGP. No other explanation 
of this behavior is known to us. An alternate explanation would have to 
address both the magnitude of the effect and its energy dependence. The 
available data is consistent with QGP; however, good experimental results 
are needed to confirm this intriguing trend. 

On the other hand, we note that the absolute yields A, and p are not 
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Fig. 2. Top: observed mid-rapidity particle yield ratio A/p as function of nucleon- 
nucleon reaction energy ^/snn- Bottom: statistical QGP hadronization total yield 
ratios in different QGP breakup scenarios. NA49 compilation of own, AGS data 
and theoretical results. 



large. At the top AGS energy (11.6 A GeV for central Au-Au collisions), 
in the most central reactions, one A or p will be produced in one out of 100 
Au-Au collisions. There are important kinetic rescattering processes which 
can generate a shift A «-> p inducing suppression, or enhancement of A or p. 
Qualitative kinetic model rescattering arguments favor A over p yield: the 
annihilation cross section on baryon matter is smaller for A. The strangeness 
exchange cross reaction is exothermic for p+ (qs) — > A+ (qq). Thus, kinetic 
rescattering in HG could increase the ratio A/p beyond relative chemical 
equilibrium yield. 
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Such dynamical effects are not common but must be investigated before 
we conclude that QGP is formed at top AGS energies When good data 
at AGS reaction energy becomes available, A/p is one of the observables 
meriting very careful experimental and theoretical study. This is the low 
reaction energy probe of QGP most resembling Jfl? suppression: the signal 
is very clean, but must be carefully considered for competing effects. 

3. Strange (anti)baryon enhancement at SPS and RHIC 

CERN experiments WA97 and NA57 have focused on the study of the 
systematics of the strange (anti)baryon enhancement with reaction energy, 
and centrality in Pb-Pb collisions. A, H and $7 and antiparticle yields have 
been measured at central rapidity y and medium transverse momentum p± 
as functions of the centrality of the collision. Comparing the yields in Pb-Pb 
to those in pBe interactions, considerable enhancement of yield per partic- 
ipant is observed. This enhancement increases linearly with the centrality, 
and geometrically with the strangeness content in hyperons, reaching a fac- 
tor of about 20 for the O + O in the central Pb-Pb collisions. Recently, 
the final results for the centrality dependence of (anti)hyperon production 
in Pb-Pb, pPb and pBe collisions at 158 A GeV/c have been reported |Hj. 
The open symbols, in figure |B1 show these results, which follow the pattern 
predicted in the recombination-hadronization model The enhancement 
rises, both with the strangeness content in the hadron, and with the par- 
ticipant number A (centrality), that is the size of the reaction region. The 
magnitude of the the enhancement is nearly the same as seen at much lower 
SPS energy range. 

There are further features of the WA97 and NA57 experimental results 
which are particularly important: the enhancement seen at 40 A GeV/c 
is very similar, nearly identical, particle by particle, as seen at 158 A 
GeV/c |101 111) expert for the most peripheral class of events. We read 
this result to mean that the peripheral 40 A GeV/c reactions do not always 
reach deconfinement. The final experimental results differ little from the 
preliminary conference reports which we discussed at this school in 2003 |12j 
and in 1999 |13| . For this reason, we focus discussion of this enhancement 
pattern at RHIC, and the comparison of RHIC with SPS. 

The solid symbols, in figure |SJ correspond to the most recent STAR re- 
sults obtained at ^/snn = 200 GeV [Tl]. The yield of (multi) strange baryons 
A(uds), 3~(dss) and antibaryons per participant N paI t in the reaction, di- 
vided by a reference yield obtained in pp reactions is shown. Within error, 
considering also the base yield, the enhancement at RHIC and SPS appears 
the same. However, there are two issues to consider: 

a) The enhancement computed by NA57 is based on pBe, where some en- 
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Fig. 3. Yields per participant AT part , for NA57-SPS Pb-Pb 
(open symbols) relative to pBe and for STAR-RHIC Au-Au 
(filled symbols) relative to pp. On left baryons, on right antibaryons and Q + O 
(triangles), circles are A and A, squares are S - and S + . Error bars represent those 
from the heavy-ion measurement. Ranges for pp and pBe reference data indicate 
the statistical and systematic uncertainty. After |14j . 



hancement of A is expected to be present, as compared to pp results. The 
NA49 SPS experiment has just reported A enhancement evaluated with 
reference to pp |15| . In this more directly comparable case, the SPS A- 
enhancement comes to be factor 5-6 (rather than 4-5) in most central Pb-Pb 
reactions. 

b) The behavior of A enhancement at SPS breaks the ranks in that it is 
seen to be nearly flat as function of centrality, and it is smaller at SPS than 
at RHIC. _ 
We conclude that the enhancement pattern of A and A is influenced by 
the prevailing baryon density, for A the enhancement is greater at RHIC 
where the baryon density is smaller than at SPS, while the reverse is true 
for A. This baryochemical potential effect is mostly erased considering d 
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S~, H , but is visible within the error bar. The enhancement of + ft is 
largest, since production of these particles is very difficult in the elementary 
reactions, especially so at the lower SPS energy. 

The interpretation we pursue is that the strange antibaryon enhance- 
ment is due to an increased yield density of strange quarks at hadronization, 
growing within the geometric source size, the increase driven by the longer 
lifespan. The gradual increase of the enhancement over the range of -/Vp art 
is an important indicator of the kinetic strangeness production mechanism. 
We return to this in section |SJ in the strangeness yield analysis as function of 
centrality |16j . This enhancement at RHIC is in fact fully understood. We 
compare the theoretical predictions with the experimental results for A and 
H~ on left and A and H + at RHIC in figure HJ Note that the normalization 
of the enhancement is different ^7], as compared to figure |21 This has oc- 
curred since the interpretation of the reference pp data changed. However, 
this vertical shift (in a log-figure) does not affect the comparison of theory 
and experiment since both use the same experimental pp yield normalizer 
to define the enhancement. We note that data in figure El is "newer". 

The theoretical lines, in figure^! follow from study of strangeness yield in 
these reactions based mainly on PHENIX single strange hadron results ^S] , 
and application of statistical hadronization model using these results |19j . 
to predict the yields shown For most peripheral collisions, this study 
is not entirely applicable, since the SHM fit to the data is not fully reliable. 
This is the 'kink' in the theoretical curves seen in figure^ An important 
help in understanding these results is the gradual growth of the enhance- 
ment with centrality. As we see, the trend is quantitatively consistent with 
strangeness yield growth (see section |f)J). Both the energy and centrality de- 
pendence of the enhancement contradicts the canonical suppression model, 
which reinterprets this enhancement as being due to reference data suppres- 
sion |2D|. 

4. From Baryons to Mesons 

The enhancement of (strange) (anti)baryons we just presented is in part 
due to enhanced production of baryons as compared to mesons when AA 
reactions are compared to pp reactions. The figure |S] shows the ratios of 
p/ir and A/Ks from central Au-Au collisions at y / sjv]v = 130, and 200 GeV 
measured by PHENIX [H] and STAR |2T1 122 1 123] . 

We recognize that the formation of mesons and baryons occurs by a dif- 
ferent mechanism than in the elementary pp reactions. In particular, the 
large baryon to meson ratio, seen in figure EJ cannot be accommodated by 
the conventional string fragmentation scheme developed for the elementary 
e + e~ and pp reactions. The quark recombination models |241 125] provided 
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Fig. 4. Yields per participant AT par t relative to pp of A and S on left and A and 
S + on right in Au+Au collisions at y / SNN = 200 GeV ^3 (compare figure O}- Our 
theoretical results |16j are for chemical non-equilibrium and are normalized to the 
same pp experimental results. 



a satisfactory description of the particle yields, in particular, the large pro- 
duction of anti baryons in the intermediate p± region. 

Within the statistical hadronization model, when chemical equilibrium 
is assumed, the relative ratio of baryons to mesons is fixed, and varies along 
with hadronization temperature T. When chemical nonequilibrium ideas 
are introduced, the parameter controlling the relative abundance of baryons 
with respect to mesons is j q , 

Baryons >y s m 

ls _ = 79 fl P (r 1 -,-), (6) 

where \ii are chemical potentials, and ^s/lq is the relative strangeness to 
light quark phase space occupancy. There is excess baryon-antibaryon pair 
yield over chemical equilibrium for j q > 1, expected when the source of par- 
ticles is a very dense deconfined quark state. Conversely, for low energies or 
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Fig. 5. Ratios of A to Ks from Au-Au and pp collisions (STAR) and p to ir from Au- 
Au collisions (PHENIX) as a function of transverse momentum In addition 

to resonance contributions in all hadrons, the A includes contributions from S° 
decays. 



small systems, we expect j q < 1. The value j q = 1 can be established when 
there is time to scatter and equilibrate the yields of mesons and baryons, 
which in general is not the case in all reactions described here. 

Consideration of chemical nonequilibrium and the possibility that in 
some reactions j q > 1, and in others j q < 1, provides an opportunity 
to describe the energy dependence of particle production, studied by the 
experiment NA49 This has as special objective the interpretation 

of the K + /7r + ratio |271 128j . which shows a pronounced energy dependent 
structure, see figure El 

The rapid rise in the K + /tt + oc su/du ratio is due to the more rapid 
increase in s than d when the baryon density is very high. The decrease 
beyond the peak is driven by a reduction of baryon stopping with increasing 
reaction energy, at which time there is a faster increase in the d yield. In 
the SHM description the key ingredient needed to obtain this behavior, in 
consistency with all particle yields, is the switch-over at the peak of the 
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Fig. 6. K + /tt + ratio of NA49 experiment [27J EE| Solid line: Result of chemical 
noncquilibrium hadron production analysis |26| . 

ratio, from j q < 1 to 7 g > 1. The gradual rise of K + /ir + at high energy is 
associated with an increase of ^s/lq- 

The observable K + /-7r + has been long considered as a signature of strang- 
eness enhancement It turned out to also be a signature of anti-quark 
suppression in high baryon density fireballs created at low reaction energy. 
Moreover, there appears to be a change in hadronization mechanism as- 
sociated with the intermediate peak in K + /tt + . The advancement of this 
most interesting physics result is mainly due to the diligent work of Marek 
Gazdzicki [27]. 



5. Strangeness per Entropy 

Individual particle yields are in general probing conditions prevailing 
at the time of QGP hadronization. However, bulk fireball properties reach 
more deeply into the history of QGP evolution. Of particular interest is the 
ratio strangeness per entropy s/S. Studies of QGP suggest that entropy 
and to a lesser extent strangeness yields are produced in the earliest stages 
of the heavy ion reaction. In the next section H3 we will show how 'deep' 
strangeness probes. In this section, we look at the qualitative features of 
this observable. 

The deconfined state is rich in both in entropy and strangeness. The 
enhancement of entropy S arises because the color bonds are broken and 
gluons can be created. Enhancement of strangeness s arises because the 
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mass threshold for strangeness excitation is considerably lower in QGP than 
in hadron matter. While entropy S determines the hadron multiplicity 
content, strangeness s determines strange hadron content. Thus, to measure 
s/S, we need to have measured particle yields which comprise much of the 
total strange and non-strange hadron multiplicity. 

In QGP, strangeness s originates primarily in gluon based reactions, e.g., 
GG — ► ss, and higher order gluon related processes. Gluons are the source 
of strangeness; their presence requires deconfinement, and is a result of 
entropy production, i.e., excess entropy S content. It is generally believed 
that the entropy production is completed after first 0.5-1 fm/c reaction 
time, at which time the parton matter has turned into the thermalized 
QGP phase. Strangeness is produced shortly after. We have been able 
to unravel much of mystery behind strangeness production (see section EJ). 
However, mechanisms of entropy formation are not understood. 

We estimate the expected magnitude of s/S. For thermally equilibrated 
non-interacting QGP phase at a temperature T s at which formation of 
strangeness essentially has ended: 

s Ps = (3/TT 2 )T3(m s /T s ) 2 K 2 (m s /T s ) = 0.03 

S ~ a (32vr 2 /45)T3 + n f [(7tt 2 /15)T3 + n 2 T s ) 1 + 0.054(ln X q ) 2 ' 

Here, we used for the number of flavors ri{ = 2.5 and 2m s /T s = 1. We see 
that the result is a slowly changing function of X q ; for large X q ~ 4, present 
at modest SPS energies, the value of s/S is reduced by 10% as compared 
to X q ~ 1 applicable to LHC and RHIC. Considering the slow dependence 
on x = trie /T ~ 1 of W(x) = x 2 K 2 {x) there is a minor residual dependence 
on the value of T s . 

The dependence on the degree of chemical equilibration which dominates 
is easily obtained by separating the different degrees of freedom: 

s 7s QGP 

8 = °-° 3 0.38 7 g GP + 0.12 7 ? GP + 0.5 79 QGP + 0.054 7g QGP (ln X q ) 2 ' ^ 

Given Eq. (jHJ), we expect to see a gradual increase in s/S as the QGP source 
of particles approaches chemical equilibrium with increasing collision energy 
and/or increasing volume. 

How does this simple prediction compare to experiment? A precise mea- 
surement of s/S is not easy, since it requires measurement of several particles 
species to considerable precision. Analysis of existent data was done as a 
function of centrality 19 , and as a function of reaction energy [^Sj. The 
centrality study included 11 centrality bins of Au-Au at ^Jsnn = 200 GeV 
for which dN/dy for 7r ± , K ± , p and p at ycM = have been presented in 
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table VIII, in Ref. ^Hj (PHENIX). These 6 particle yields and their ratios 
change rapidly as function of centrality. In addition results from STAR for 
K*(892)/K~ |30| . and </>/K~ which show minor centrality dependence 
were used. 

Supplemental constraints (ratio tt + /ir~ = 1. db 0.02, strangeness conser- 
vation, electrical charge to net baryon ratio) help to determine the best fit. 
The 7 SHM parameters (volume per unit of rapidity dV/dy, temperature 
T, four chemical parameters A g , A s , 7 g , 7^ and the isospin factor A/3 are in 
this case studied in a systematic fashion as function of impact parameter 
using 11 yields and/or ratios and/or constraints, containing one (pion ratio) 
redundancy [T5] . 

X 2 minimization yielding good significance is easily accomplished, show- 
ing good consistency of the data sample. The resulting statistical param- 
eters are seen in figured on left, as function of participant number. The 
full non-equilibrium model 7 9 / 1,7 S / 1 (full circles, blue) and semi- 
equilibrium j„ = 1 (open circles, red) are shown. From top to bottom 
the (chemical) freeze-out temperature T, the occupancy factors 7 g , 7^/7^ 
and together in the bottom panel the baryon hb an d hyperon fig chemical 
potentials. 

Each of the panels shows very interesting behavior indicating some par- 
ticular physics feature at chemical freeze-out. Perhaps the overall most 
striking result is that, for A > 20 only, 7^/7^ varies significantly, as does 
(not shown here) the volume of the reaction domain dN/dV. The scaling 
of dN/dV with A 2 /3 (TBI is expected in a geometric reaction picture. The 
growth of 7s /7q indicates that with size of the system there is more effec- 
tive cooking of strangeness. There is no saturation of j s as we approach 
the most central reactions. This is inherent in the data we consider which 
includes the yields of (f>, and K*. This result is consistent with the impli- 
cation that strangeness is not fully saturated in the QGP source, though it 
appears over-saturated when measured in the hadron phase space. 

The baryo-chemical potential and hyperon potential (bottom panel) is 
model independent and its nearly constant value indicates that only the 
reaction energy and not the centrality of the reaction establishes baryon 
density at hadronization. We find a hadronization temperature which, when 
we allow for j q / 1, is at the level of T = 140 MeV, and it rises to T ~ 155 
when we force light quark chemical equilibrium on the particle yield fit. 

On the right, in figure [7J we show the same statistical variables as a 
function of reaction energy. To study the energy dependence, we must 
assemble several different experimental results from different facilities and 
experiments [55]. The lowest energy result is from our AGS study jHlj; the 
SPS data we used are from NA49 energy dependence exploration at the 
CERN-SPS JUj. These results are for the total particle yields. The two 
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Fig. 7. From top to bottom: temperature T, light quark phase space occupancy j q , 
the ratio of strange to light quark phase space occupancies J s /jq an d the chemical 
potentials (B for baryochemical fis and S for strangeness (is) The lines guide the 
eye. Left: as function of inelastic (wounded) participant number A. The full non- 
equilibrium model 7 9 ^ 1,7 S 7^ 1 (full circles, blue) and semi-equilibrium y q = 1 
(open circles, red) are shown. Right: energy -y/sNN-dcpcndcnce in full chemical 
nonequilibrium analysis. Note that the additional results at RHIC (two highest 
energy results in green on right) apply to central rapidity region dN/dy study, 
while the AGS/SPS and Brahms- RHIC results are for total yields (blue lines). 



highest energy results are based on studies of RHIC data at 130 and 200 
GeV at central rapidity and address the dN/dy particle yields, with the 
highest point corresponding to the results presented at greatest centrality 
on the left-hand side in figure 

There are several energy dependent features seen on the right-hand side 
of figure[7| We note the dotted vertical line at energy where some of the SHM 
parameters undergo a change of behavior. The value of energy corresponds 
to where NA49 reports a change in the behavior of K + /-7r + oc s/d. 

The statistical parameters obtained in the analysis are used to compute 
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s/S as function of centrality ^H]) and as function of reaction energy |2(i| . 
One indeed finds that corresponding to the growth of Js/lq seen in figure [7[ 
the value of s/S grows, reaching nearly the QGP chemical equilibrium value, 
s/S ~ 0.03 at highest RHIC energy in most central reactions. As a function 
of centrality at top RHIC energy flat peripheral behavior where s/S ?S 0.02 
at A < 15 turns into smoothly increasing s/S reaching s/S ~ 0.03 in most 
central reactions. As a function of energy near to -^/snn — 6.5 GeV in most 
central collisions, there is a change in growth rate of s/S. 

The significance of s/S as an observable of QGP is that it is not depen- 
dent on the details of the dynamics of hadronization which remain unre- 
solved today. The value is predictable, and helps us understand the degree 
of chemical equilibration of QGP prior to hadronization. The numeric value 
for saturated QGP phase is directly proportional to the ratio of the num- 
ber of strange degrees of freedom to all degrees of freedom in the plasma 
(s/S oc g s /4g, see next section). Thus, s/S ~ 0.03 indicates that quark- 
gluon quanta are active. The smooth rise and the magnitude of s/S are 
consistent with the QGP reaction picture — a full interpretation of s/S 
requires better understanding of QGP chemical equilibration, a topic we 
address next. 



6. Strangeness production at RHIC and LHC 

The kinetic description of strangeness production can be reorganized 
to directly evaluate the strangeness yield normalized by entropy yield in 
plasma, s(T(t)/ S(T(t)). In this way, the model dependence on the ex- 
pansion-dilution phenomenon due to rapid expansion of the QGP phase is 
considerably reduced. We describe our recent progress in the study of the 
chemical equilibration at RHIC and LHC [33], and evaluate how this work 
supports the interpretation of the results of data analysis for s/S presented 
above. 

The temporal evolution of s/S, in an expanding plasma, is governed by: 

When all ji — > 1, the Boltzmann collision term vanishes, the system has 
reached chemical equilibrium. The value arrived at for the observable s/S 
depends on the history of how the system evolves at early stages, and reaches 
QGP chemical equilibrium for gluons in particular. 

In order to be able to solve Eq. Q, a relation between s/S and 7 S is 
required: 

s 9s 90 2 / ka s (T) 



8=^-8^* K *(*)> Z = m s/T, g s = 2 s 3 c ll ^ ■ ■ ■ ) ■ ( 10 ) 
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In the initial period, gluons and quarks have not reached chemical equilib- 
rium, thus the actual numerical integrals of Bose and Fermi distributions 
are dependent on the values 7 9ig . The quark-gluon QCD interactions are 
introduced to describe the effective degeneracy with precision: 

50a s (T) \ 
1 — . . . J + strangeness. 

2l7T J 

(11) 

The strangeness component is made consistent with g s in Eq. ()10j) and the 
concurrently numerically computed 7 S . 

The degeneracies, seen in Eq. (jlOjl for the entropy, include the effect of 
interactions. Therefore, consistency requires that the interaction effect is 
also introduced in the strange quark degeneracy g s , Eq. (|10[). The value of 
k = 2 applies to massless strange quarks. At T = (or said differently, for 
m^> T) the early study of quark matter self-energy suggests that k — ► |34| . 
The full range < k < 2 has been explored, the results for k = 1 are 
presented. 

We need to model the expansion dynamics of QGP matter. The study of 
hadronization at RHIC provides a strong constraint on the initial condition 
and the evolution dynamics of QGP at RHIC. For LHC three modifications 
are introduced: 

1. To account for the greater reaction energy, we increase the initial 
entropy dS(T )/dy by factor 4 from 5000 at RHIC to 20,000 at LHC. We 
do not change in the initial value of s/S\ TQ = 0.016; thus, in elementary 
parton-string interactions, the relative strength of strangeness and non- 
strange hadron production is left unchanged. This implies an increase in 
initial strangeness yield by a factor 4 at LHC compared to RHIC. 

2. The initial parton thermalization time is reduced from To = 1/4 fm 
at RHIC to To = 1/10 fm at LHC. This does not influence the outcome. 

3. In order to accommodate the greater transverse expansion pressure, 
we increased the maximum transverse flow velocity reached to v± = 0.80c . 
A higher value would increase the over saturation of QGP phase space. 
Despite a much greater expansion velocity, the evolution time at LHC is 
significantly longer compared to LHC, with the most central collisions tak- 
ing up to 30% longer to reach the freeze-out hadronization temperature 
Tf = 0.14-0.17 GeV. This is due to the greater initial entropy and energy 
content that takes longer to dilute to the hadronization condition. 

Finally, we need to evaluate the invariant strangeness production rate 
per unit 4- volume, A"^ ss and A q ^ ss , see Eq. @. These are obtained 
employing the running strength of the QCD coupling, with a s (/i = m^o) = 
0.118 evolved to the applicable energy domain \x using the two loops (3- 
function. 



g = 2 s 8 c 7 g 1 



15a, (T) 

47T 



+-2 s 3 c 2f7g 
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Strangeness production, a relatively soft process, can be described per- 
turbatively for two reasons: 

a) the reaction processes which change yield of strangeness can compete with 
the fast v± > 0.5c expansion of QGP only for T > T s = 220 MeV ~ 2m s , 
for lower temperatures the strange quark yields effectively do not change. 
T s is the strange quark chemical freeze-out temperature in QGP. Using the 
relation fi = 2irT, this implies that all strangeness yield evolution occurs 
for n > 1.4 GeV. 

b) Because of the magnitude a s (fi = m z o) = 0.118, one can run a s to the 
scale of interest, [i > 1.2GeV, and the value cx s /tt £ 0.2 allows perturbative 
QCD methods. 

We also evolve the central strange quark mass value m s (n = 2 GeV) = 
100 ± 25 MeV. We use as the energy scale the CM-reaction energy [i ~ \fs 
working in two loops. Some simplification is achieved by considering at 
temperature T the QCD scale [i ~ 2ttT, thus m s (T) = m s ([i = 2ttT) with 
m s (T = 318 MeV) = 0.1 GeV. The actual temperature T(t), and thus the 
time dependent values of the strange quark mass, are seen in the top panel 
of figure 03 

The figure El illustrates in detail how a variation of the unknown initial 
chemical equilibration conditions impacts the final strangeness yield results. 
Since the entropy seen at the end of the reaction is already present at an 
early time, it is fixed as function of time and is the same for the different 
initial state scenarios explored. T(to) is obtained (top panel) for a wide 
range of assumed initial gluon (and quark) occupancy 0.1 < 7 g (To) < 2.1, 
in step of 0.5, shown in the middle panel by dashed lines. 

Results for strangeness production at RHIC (left) and LHC (right) as 
function of proper time r are shown in middle and bottom panels; the solid 
lines in middle panel show the resulting 7 S , and in the bottom panel s/S. 
Two results emerge: the value of 7 S at QGP hadronization is nearly inde- 
pendent of the initial conditions; the evolution of s/S in time is a universal 
curve. A large part of the difference between RHIC and LHC arises because 
the expansion lasts longer. At RHIC, the thermal strangeness production 
raises the value of s/S from 0.016 to 0.028, and at the LHC the thermal 
strangeness production raises s/S from 0.016 to 0.032. This 20% increase 
in relative strangeness yield has a highly significant impact on production 
of strange hadrons. 

These results show that the observable s/S is dependent on the dynam- 
ics of dense matter evolution (and also QCD parameters a s ,m s , but not 
on the internal structure of the QGP, and in particular not on the degree 
of chemical equilibration). At a fixed entropy yield, a plasma that has not 
yet reached chemical equilibrium is hotter. This higher T compensates in 
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Fig. 8. Strangeness production at RHIC (left) and LHC (right) for the case of 
5% most central collisions. s/S (bottom panel) and 7 S (solid lines middle panel) 
evolution as function of proper time r. In this study, widely varying initial gluon 
conditions (7 g dashed middle panel) constrained to same entropy content, and thus, 
varying temperature T (top panel), is presented. 

strangeness production the lack of chemical equilibration, i.e., 7g < 1. En- 
tropy is practically constant since the process of chemical equilibration does 
not produce significant amount of entropy jSS]) which is mostly produced 
in the process of thermal equilibration. 

The final strangeness content expressed by s/S depends on the initial 
entropy content, dS(To)/dy. This value can change both due to change in 
energy and centrality of the reaction. A scenario study shows that variation 
of centrality leaves us on the universal s/S curve seen in figure 03 but the 
freeze-out temperature is reached sooner for smaller volumes with smaller 
dS/dy, and thus, smaller s/S is achieved. A similar result is obtained when 
varying reaction energy. For this reason, we present in figure|§]the final s/S 
and 7 S as function of dS/dy. These two results, s/S and j g , will change 
independently under variation of strange quark mass, with 7 S remaining 
constant and s/S changing. Considering extreme conditions for the volume 
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Fig. 9. s/S (top panel) and j s (bottom panel) as function of dS/dy. 



evolution one can vary both s/S and 7 S by the same magnitude, 10-20%. 
We remind the reader that the 7 S discussed in this section is for the QGP 
phase while the 7 S discussed in hadronization analysis is for HG phase and 
is generally 2.5 times greater. 



7. Strange Hadron Resonances 

Theoretical models of hadron production need to account in their study 
of particle spectra or/and hadron yields for the contribution originating in 
resonance decays. In the statistical hadronization, the production of hadron 
resonances is abundant. Within a particle 'family', particle yields with same 
valance quark content are in relation to each other well described by inte- 
grals of relativistic phase space, eventually corrected for rescattering |36| . 
The relative yields of, e.g., K*(sq) and K(sq), are controlled by the particle 
masses mj, statistical weights (degeneracy) gi and the hadronization tem- 
perature T. In the Boltzmann limit, one has (star denotes the resonance): 



N* 
~N 



g*m* 2 K 2 (m*/T) \M*\ 
gm 2 K 2 (m/T) \M\< 



(12) 



In the statistical hadronization model, the assumption made is that the 
matrix element A4(E) is not strongly energy dependent (E is energy of 
particle produced). Thus, as shown above, the energy average factorizes. 
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Moreover, one also expects that for most strongly interacting particles, the 
quantum strength saturates at its maximum, \Ai\ 2 — > 1, |A^*| 2 — ► 1, in the 
energy range in which most particles are produced (but not necessarily in 
some extreme corners of phase space where special situations prevail). 

Resonance production is measured, in an experiment, directly by the 
decay product invariant mass method. Presence of these particle, and the 
fact that their 4- momenta add up to form invariant mass of the parent means 
that most if not all of the decay products do not significantly rescatter. 
Except for special cases ((f), f2, U) resonances contribute 50% and more of the 
yield of each hadron. When this is recognized in an analysis of experimental 
data, one finds that all particles are produced at the same temperature and 
that the yield temperature (chemical freeze-out) and spectral temperature 
(thermal freeze-out) are the same. This analysis was done for the maximal 
SPS energy |37| . and the maximal RHIC energy jHH]. Claims to the contrary 
rely on analysis which rescatter the resonance decay products such that all 
particles considered are forced to satisfy a thermal spectrum overlaid with 
flow from expansion dynamics. This, then, produces for different particles 
different freeze-out conditions of temperature and velocity. Such analysis is 
of course wrong as it presumes tacitly that not a single resonance should be 
observable by invariant mass method. 

Putting aside the global analysis of data, we can argue qualitatively 
for single freeze-out, comparing the baryon and antibaryon mj_ spectra as 
measured very precisely by WA97 and NA57 at SPS [SHI EH], and STAR at 
RHIC |411 142j . We note that if hadron gas were to surround at any relevant 
stage the plasma phase, such a symmetry could not arise in the presence 
of a baryon excess: the A, 3 annihilation is strongly momentum dependent, 
this should deplete the low momentum yield. This deforms the shape of 
the antihyperon spectra, as compared to the spectra of hyperons. More 
specifically, 3 could annihilate on relatively abundant N — N, A — A, X — X 
(rather than the much smaller 3 — 3), and yet the spectral mj_-distributions 
of 3 are within a fraction of percent error, identical to 3. Absence of 
spectral deformation suggests that 3 emerge directly and without significant 
rescattering. 

Our working hypothesis is therefore that hadronization of the QGP de- 
confined phase formed in high energy nuclear collisions at top SPS energy 
and RHIC energies is direct, fast (sudden) and occurs without significant 
sequel interactions. In this circumstance, the measurement of the relative 
yield of hadron resonances is a sensitive test of the statistical hadronization 
model (SHM) hypothesis and lays the foundation to the application of the 
SHM method in data analysis. One of the great successes of the chemical 
nonequilibrium particle yield analysis is the agreement with the experimen- 
tal resonance yields. These, in particular, are at RHIC K*(892), A(1520) jlSl 
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144] . Our relatively low hadronization (chemical freeze-out) temperature 
makes this possible. More data is on the way and will keep us busy in 
coming years [451 146j . 

We believe that the systematic study of resonance production can re- 
solve the question about the reaction energy dependence of the freeze-out 
condition temperature. One could argue naively that this is simply the 
temperature at the quark-hadron transformation seen in lattice computa- 
tions. This implies a smooth change with reaction energy, since temperature 
rises smoothly while baryon density decreases. The maximum value of T is 
reached at highest reaction energy. However, since there is no sharp phase 
transition at small baryon density prevailing at LHC, RHIC and perhaps 
top SPS energy, there is in this hypothesis a lot of uncertainty: the value of 
cross-over temperature varies depending on observables considered. More- 
over, among lattice results there is considerable difference in the range of 
interest for maximum T, with T = 195 ± 10 MeV jlTj. and T = 151 ± 5 
MeV [48] ■ at zero baryon density. In recent years, T has varied within these 

limits mi EE eu • 

Our view is that freeze-out temperature is, to a large extent, influenced 
by the dynamics of the rapidly expanding fireball 52_ : a) When the collec- 
tive flow occurs at parton level, the color charge flow, like a wind, pushes out 
the vacuum, adding to thermal pressure a dynamical component [531 154] . 
This can, in general, lead to supercooling and a sudden breakup of the fire- 
ball. We find that this can reduce the effective hadronization temperature 
by 20 MeV 

b) We also must consider the chemical composition of the hadronizing QGP. 
It is well known that the critical temperature is dependent on the flavor 
content. Thus, the degree of chemical equilibration of strange quark fla- 
vor influences the exact location of the boundary, and the nature of the 
transition. Because the degree of equilibration in the QGP depends on the 
collision energy, as does the collective expansion velocity, we cannot expect 
a simple hadronization scheme appropriate for the hadronization of nearly 
adiabatically expanding Universe. We can be nearly sure that the chemical 
conditions matter and can displace the transition temperature. 

The question of considerable importance is how we can find agreement 
on what the temperature of hadronization is, and be able to trace in the 
phase diagram the hadronization boundary, in a way that leads to general 
agreement. The disagreement is considerable. In figure ITU1 we see the situ- 
ation proposed for the hadron system considered to be in absolute chemical 
equilibrium 

At highest heavy ion reaction energy, one obtains chemical freeze-out 
temperature T ~ 170 MeV. Values as low as T ~ 50 MeV are reported 
at lowest reaction energies available. As the collision energy increases, the 
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Fig. 10. Freeze-out boundary temperature T and baryo-chemical potential /ib on 
reaction energy following from the chemical equilibrium hypothesis, after Ref. |55|. 
The line guides the eye. 

freeze-out temperature increases and the baryonic density (here baryonic 
chemical potential //b) decreases. Such an increase of freeze-out temperature 
with collision energy is expected on general grounds, since with increasing 
reaction energy a greater fraction of the energy is carried by mesons created 
in the collision, rather than pre-existing baryons |56j . 

There is no sign of any structure, in figure ITU1 when the reaction energy 
varies. The rather smooth hadronization curve makes it impossible to de- 
scribe the experimental discontinuity in the ratio K + /tt + (see figure EJ) and 
other related experimental features (27]. An effort was made to interpret 
this in terms of a shift from baryon to meson dominance [H7j of the hadron 
yields. However, in the chemical equilibrium model, all observables includ- 
ing K + /7r + remain a smooth function of reaction energy, in contrast to the 
experimental results. 

The systematic behavior of freeze-out condition in the T-/xb plane al- 
lowing for j q ^ 1 is quite different , see figure 1111 The two higher T 
values at right are for 20 A GeV (lowest CERN-SPS energy) and (furthest 
to right) 11.6 A GeV (highest BNL-AGS energy) reactions. In these two 
cases the source of particles is relatively large, dilute and strongly chemi- 
cally under-saturated, and hotter than otherwise has been expected. Such 
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Fig. 11. Dependence of freeze-out temperature T and baryo-chemical potential /ie 
on reaction energy in the non-equilibrium case [55]. The dashed lines guide the 
eye. 

a system could be a conventional hadron gas fireball that had not the time 
to chemically equilibrate. Other options were considered such as a phase of 
constituent massive quarks. 

The higher temperature and lower density imply that at the time of 
freeze-out the fireball is not expanding as fast as found at higher reaction 
energies. Thus without doubt sudden hadronization does not apply here. 
The transition to the supercooled regime occurs and the phase space changes 
from chemically under-saturated (j q < 1) to chemically over-saturated (j q > 
1) where the two dashed lines meet in figure 111! At heavy ion reaction 
energy below (i.e., to right in figure lllj) of this point, hadrons freeze-out 
with meson number enhanced with reference to antibaryon yield. Above 30 
A GeV the antibaryon to meson yield is enhanced. 

The relative resonance yield Eq. I|12|) is only dependent on the tempera- 
ture of hadronization. Hence the relative yield of resonances will distinguish 
the two scenarios presented in figures ITU1 and ITTI 58 . Torrieri et al. evalu- 
ated the relative yield of resonances along the two hadronization curves as 
function of reaction energy and have found that the behavior is practically 
opposite. Relative resonance yields in the equilibrium model rise monoton- 
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ically with energy of heavy ion reaction. The nonequilibrium model implies 
overall a reduction of T with increasing heavy ion reaction energy, with a 
strong dip near to the 'horn' energy. Moreover, and beyond the argument 
presented in |58| . the visible resonance yield could be negligible at lowest 
SPS energies due to rescattering of decay products in absence of sudden 
hadronization. 

8. Outlook 

Rereading many experimental and theoretical works we have had im- 
pression that there is more experimental data than theoretical effort in this 
very rich field. To our despair, much of the data analysis is scope-restricted 
and does not seek overall consistency. Authors of failed theoretical ideas 
do not come forward to withdraw. Their ideas live on, and are resurrected. 
Today, we see more confusion in the field than there has ever been. Unfor- 
tunately, this confusion is welcome in some quarters, simply because there is 
such good strangeness data, generating a very clear message, overshadowing 
other observables. 

At SPS and RHIC, we see a chemically equilibrated strangeness rich s-q- 
g-system. What we are not 100% sure of is, if this is every-day quark-gluon 
plasma, though there is plenty of evidence that it is. We have discussed how 
the degree of chemical equilibration in QGP varies, with system size, and 
reaction energy. The result follows from glue-based strangeness production. 
Entropy content seen in particle multiplicity, is for the most central and most 
energetic RHIC reactions just as QGP predicts — the ratio s/S measures 
directly the ratio of strange to all QGP degrees of freedom. At the 10% level, 
or better, it agrees with what one expects from the chemically equilibrated, 
interacting quark-gluon fluid. 

Hadronization of the fireball of matter formed in heavy ion reactions 
leads to quite different spectra and yields of hadrons than we expect based 
on elementary pp reactions. This change in reaction mechanism favors, in 
particular, the production of multi strange baryons and antibaryons. The 
enhancement we see is what statistical recombination of quarks predicts, 
both as function of centrality and energy. 

Many years ago, when the first ideas how to look for quark-matter were 
born, it seemed that the energy threshold to deconfinement could be low. 
Nucleons in nuclei seemed to be just barely retaining their individual iden- 
tities. Just a little push and squeeze with a few GeV beam could perhaps 
be sufficient to lead to deconfinement 0J. This is where we are today, after 
a long and dramatic excursion to the very high RHIC energies: Strange 
antibaryon enhancement suggests that at least down to 40 A GeV we have 
s-g-g-matter. The simplest of all possible observables, the K + /ir + ratio 
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shows a threshold between 20 and 30 A GeV projectile energy. 

Let us cross our fingers that RHIC can and will run at a few GeV per 
beam. The existent detectors would in this environment produce very pre- 
cise and have nearly full coverage in phase space for the data on strange 
hadron production, including resonances. This would with certainty resolve 
any doubt about QGP, both its formation and threshold as function of cen- 
trality and reaction energy. This will further lead to detailed understanding 
of the phases of QCD. This work will complement the LHC based study of 
perturbative QGP. 
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